OLD ISOLATED ACCRETING NEUTRON STARS: 
CONTRIBUTION TO THE SOFT X-RAY BACKGROUND 

IN THE 0.5-2 KEV BAND 



Silvia Zane , Roberto Turolla 2 , Luca Zampieri 1 , Monica Colpi 3 and Aldo Treves 1 

1 International School for Advanced Studies, Trieste 
Via Beirut 2-4, 34014 Trieste, Italy 

2 Department of Physics, University of Padova 

Via Marzolo 8, 35131 Padova, Italy 

3 Department of Physics, University of Milano 

Via Celoria 16, 20133 Milano, Italy 



1 



ABSTRACT 

The issue of the observability of Old Isolated Neutron Stars (ONSs) accreting 
from the interstellar gas, is reconsidered using the spectra presented in Zampieri 
et al. (1995). In particular, we focus our attention on the overall soft X-ray 
emission of ONSs, which may provide a substantial contribution to the X-ray 
background in the 0.5-2 keV band. It is found that if accretion is channeled 
into polar caps by a magnetic field of ~ 10 9 G, up to ~ 12 — 25% of the 
unresolved soft excess observed at high latitudes in the X-ray background can 
be explained in terms of emission from accreting isolated ONSs. However, this 
would imply the detection with ROSAT of ~ 10 sources deg -2 above a threshold 
of ~ 10~ 14 erg cm _2 s _1 . A brief reconsideration of the observability of isolated 
ONSs as individual X-ray sources is also presented. 

Subject headings: accretion, accretion disks - background radiations - stars: 
neutron - X-rays: sources - X-rays: stars 
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1. INTRODUCTION 

Old neutron stars (ONSs), i.e. neutron stars which evolved beyond the pulsar 
phase, may be as many as ~ 10 9 in the Galaxy. Even if their number accounts 
for a few percent of the total Galactic star population, isolated ONSs have not 
been unambiguously detected so far. Although the thermal emission resulting 
from cooling of the hot interior is too weak for being observed after a time 
~ 10 Gyr, old isolated neutron stars fed by the interstellar gas may show up 
as very weak, soft X-ray sources, as originally suggested by Ostriker, Rees 
and Silk (1970). For a speed relative to the interstellar medium (ISM) of 10 
kms -1 and a density of 1 cm -3 , the accretion luminosity L of an isolated ONS 
is ~ 10 31 ergs _1 . If ONSs emit from their entire surface (of radius i?*) as a 
blackbody, the effective temperature 

( L V' 4 

turns out to be ~ 30 eV. A harder spectrum would be emitted if ONSs retain, 
after ~ 10 Gyr, a residual magnetic field which funnels the accretion flow in 
the polar cap regions: for a field of 10 9 G, T e ff could be 3 times larger at the 
same accretion luminosity. The low bolometric luminosity and the softness of 
the spectrum explain the difficulty of observing an isolated accreting ONS. 

The detection of ONSs was included as a possible target for the Einstein 
mission (Helfand et al. 1980), but it was only a decade later that this issue 
came to life again when Treves and Colpi (1991, hereafter TC) reconsidered the 
observability of ONSs with ROSAT. Assuming a blackbody spectrum, using 
the velocity distribution proposed by Paczyhski (1990) and a local density 
0.07 cm -3 for the ISM, they found that thousands of ONSs should appear in 
the ROSAT PSPC All Sky Survey in the most favorable case of polar accretion. 
A complete analysis by Blaes and Madau (1993, hereafter BM) essentially 
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confirmed TC results. 

Recent observations with the ROSAT Wide Field Camera during survey 
phases (Pounds et al. 1993) and with EUVE have provided evidence for the 
presence of a relatively small number (~ 23) of unidentified bright sources: 
this led Madau & Blaes (1994) to rule out there being as many as 10 9 isolated 
ONS in the Galaxy, if accretion is spherical. Uncertainties on the statistical 
properties of the old neutron star population and on the physics of the accretion 
process make nonetheless the issue of the total number of Galactic ONSs still 
controversial (see e.g. Treves et al. 1995). The search for favorable sites where 
accretion can make a lone neutron star shining is therefore of importance. In 
connection with this problem, the detectability of ONSs in molecular cloud 
was examined by BM and Colpi, Campana, & Treves (1993, hereafter CCT). 
These studies were recently corroborated by the work of Stocke et al. (1995), 
in which it is proposed that a still unidentified source in the Einstein Extended 
Medium Sensitivity Survey is an ONS in the Cirrus cloud. A systematic search 
of isolated ONSs using ROSAT observations is presently under way (Danner, 
Kulkarni & Hasinger 1994). 

The interest for the X-ray emission from isolated ONSs is not restricted 
to the detection of individual sources. Despite their intrinsic weakness, in fact, 
they could provide a contribution to the diffuse X-ray background (XRB), being 
their total number very large. This point was already discussed by BM, who 
found, however, that the integrated flux from ONSs is much lower than the 
XRB above 0.5 keV. In a recent paper Hasinger et al. (1993) have analyzed 
the X-ray emission in selected high-latitude fields and found that the number 
of resolved sources exceeds by ~ 60% the number of QSOs expected from a 
standard evolutionary scenario. This led the suggestion that a new population 
of sources can contribute sizably to the soft X-ray background detected at 
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0.5-2 keV with ROSAT. Although such a population may be extragalactic in 
origin, Maoz & Grindlay (1995, hereafter MG) showed that its properties are 
compatible with those of galactic objects and tentatively identified them with 
Cataclysmic Variables (CVs). Accreting ONSs could have the desired spatial 
distribution but were ruled out as possible candidates by MG on the basis of 
their too soft emission. We stress, however, that both BM and MG restricted 
their discussion to the very simple case in which the spectrum is a blackbody 
at T e ff, emitted from the entire star surface. 

In the light of the previous discussion, the issue of the correct determination 
of the spectrum of an accreting ONS seems indeed to be compelling. In a recent 
paper (Zampieri, Turolla, Zane & Treves 1995), we addressed this problem 
presenting a self-consistent calculation of the X-ray spectrum emitted by a 
unmagnetized neutron star, solving the full radiative transfer problem for a 
static, plane-parallel atmosphere. It was found that, at very low accretion rates, 
corresponding to L^~> 10 31 ergs -1 , emitted spectra show an overall hardening 
with respect to the blackbody at T e //, even if the atmosphere is in LTE and 
the scattering optical depth is lower than the absorption one; typical hardening 
factors (( radiation temperature )/T e ff) are about 2.5. At energies larger than 
0.5 keV, where the blackbody emission becomes vanishingly small, the computed 
spectrum shows a persistent tail and a fraction of 10%-20% of the total flux is 
emitted above this threshold. 

In this paper we reconsider the observational impact of accreting ONSs 
using the calculated spectra and allowing for polar cap accretion. The paper is 
structured as follows: in section 2 we overview the main points about the physics 
of the accretion process and in section 3 the present distribution function of 
ONSs is calculated. Section 4 contains our analysis of the contribution of ONSs 
to the XRB and in section 5 estimates of the observability of individual sources 
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are presented. Discussion and conclusions follow. 

2. ACCRETION ONTO ISOLATED NEUTRON STARS 

The total luminosity L emitted in the accretion process depends on the star 
velocity v and the density of the surrounding gas, n, and can be written as (see 
e.g. Novikov & Thorne 1973) 

£ = (1 -y)m = 1.85 x 1(T 3 (1 - y)nv~ 3 (1) 

where n = (1 + x) n H, X = 0.36 for standard chemical composition, n# is 
the hydrogen number density, £, m are the luminosity and the accretion rate 
in Eddington units, and the velocity is in km/s. The factor 1 — y, where y = 



a/(1 — Rg/R*) and R g is the gravitational radius, is the relativistic efficiency; 
typical values for the neutron star mass and radius are used throughout: M = 
1.4M* and R* = 3R g = 12.6 km. In the previous formula the star velocity 
is replaced by the local sound speed of the interstellar gas c s when v < c s . 
Alternatively one could assume that £ oc (v 2 + c^) -3 / 2 . It should be noted, 
however, that when v ~ c s , both approaches give an approximated value, no 
simple expression for m being available (see again Novikov & Thorne 1973). 

In order to evaluate the accretion luminosity a value for nu must be 
supplied. The spatial distribution of the interstellar medium in the galaxy is 
highly inhomogeneous (see e.g. Dickey & Lockman (1990) for a detailed review 
on this subject). Observational data from Lya and 21 cm absorption measures 
shows that the ISM distribution of both cold and warm HI is nearly constant 
in radius while its ^-dependence can be fitted by 

n HI = n t exp (-^ + n 2 exp + n 3 exp (-£) (2) 

with m = 0.395, n 2 = 0.107, n 3 = 0.064; a x = 212, a 2 = 530 and h = 
403 pc; rii and Oi are in cm -3 and pc, respectively. The applicability of the 
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previous expression is restricted to the range 0.4 < R/Rq < 1, where R is the 
galactocentric radius and Rq = 8.5 kpc is the distance of the Sun from the 
Galactic center. The gas layer has a scale height of about 230 pc in the vicinity 
of the Sun, while for R < 0AR o it shrinks to « 100 pc and in the outer Galaxy 
it expands linearly up to ~ 2 kpc. The other important contribution to the 
total ISM density comes from molecular hydrogen. The best tracer of H 2 is the 
CO molecule, and observational data suggest a local gaussian distribution with 
a scale height of ~ 60 — 75 pc. Observations, however, are much less conclusive 
as far as the midplane density is concerned (Bloemen 1987), also because it 
may significantly depend on R (De Boer 1991). Here we approximate the H 2 
distribution with a gaussian with central density 0.6 cm -3 and FWHM 70 pc 
(De Boer 1991). The ionized component gives only a very small contribution 
to the total density, and it will be neglected. 

The remaining free parameter in our discussion is the value of c s . If the 
star is moving supersonically through the ISM, the main contribution to the 
accretion rate comes from the material swept by its motion. On the other 
hand, we are mainly interested in the low-velocity tail of the ONSs distribution 
because these are the stars which radiate the higher luminosities. Their velocity 
may become lower than the ISM sound speed, and the accretion rate is then 
fixed by the thermal velocity. Since the accretion radius turns out to be always 
smaller than the Stromgren radius by a factor 10 3 , c s is the sound speed in a 
fully ionized medium at the equilibrium ionization temperature (T ~ 10 4 K) 

c s ~ lOv^T/lO 4 !^) kms" 1 ~ 10 kms" 1 . (3) 

Once nn and c s are fixed, the observability of an ONS depends on its 
distance, velocity and on the response of the detector. The count rate N, 
corrected for the absorption due to the interstellar gas and for the detector 
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effective area A u , is 

f°° F 

N = J -^-exp(-a„N H )A u dv (4) 

where F v = L^/Aitd 2 , d is the distance of the source, N H is the column density 
and a v is the absorption cross section (Morrison & McCammon 1983). 

Spectral calculations are presented in Zampieri et al. (1995) for emission 
coming from the entire star surface and for 10 -7 < £ < 10~ 2 . The lower limit 
corresponds to L ~ 2 x 10 31 erg/s for M* = 1.4M and a pure H atmosphere. 
Since luminosities ~ 10 30 erg/s are expected from ONSs with velocities ~ 20—30 
km/s, we have calculated model spectra down to £ = 4 x 10 -8 ; numerical 
problems prevented us to reach lower luminosities. Spectra corresponding to 
£ = 10~ 7 , 4 x 10~ 8 are shown in figure 1. A comparison between the actual 
emerging flux and the blackbody one above 0.1 keV shows that the former is 
greater by 10%-40%, and that it deviates more and more from a Planckian 
as luminosity decreases, showing a broad maximum. In the following sections 
we will be interested to evaluate the source luminosity in the 0.5-2 keV band, 
L x , for any given star velocity v, or, equivalently, for any given value of the 
bolometric luminosity L. This is done using the analytical expression 

j = 1.98 + 0.26 log i (5) 

which fits the values of £ x for numerical models to better than 4 %. 

We account for the presence of a relic magnetic field assuming that 
emission is now concentrated in the polar caps. All radiative effects of the 
magnetic field are neglected, so its only role is to channel the accretion flow. 
A different approach was considered by Nelson et al. (1995), who evaluated 
the contribution of the cyclotron line and supposed that the thermal part of 
the emission is blackbody. For very high fields, the atmosphere is no more in 
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LTE and a detailed analysis of radiative transfer including magnetic effects is 
required (see e.g. Miller 1992). In addition, in this case magnetic pressure 
could temporarily stop the incoming flow, and the accretion process may be 
cyclic with recurrence time ~ 10 6 s (Treves, Colpi & Lipunov 1994). For this 
reason, here we limit our discussion to the case of B ~ 10 9 G. The total emitting 
area is now 2A C , where A c = irr^ is the area of the polar cap, r c = R^J 2 r^^ 2 
and ta is the Alfven radius: 

/ B 2 Rt \ 2/7 



T A = 



\V2GMMJ 

OxXOW'^G^f (g^r- (6) 



The limited size of the radiating region produces a hardening of the 
spectrum with respect to the unmagnetized case with the same luminosity, 
since now it is F = L/2A C > L/AttRI- 

3. THE NEUTRON STARS DISTRIBUTION FUNCTION 

ONSs accretion luminosity strongly depends on the star velocity with respect to 
the ISM and on the local density of the surrounding material, so, in order to shed 
light on their collective observational properties, the knowledge of the present 
ONSs distribution function in phase space is needed. Unfortunately, the velocity 
distribution of pulsars at birth is still poorly known, and any evolutionary 
scenario remains affected by this indetermination. In a detailed study, Narayan 
& Ostriker (1990) found that observational data of periods and magnetic fields 
for a sample of about 300 pulsars are well fitted by invoking the presence of 
two populations of neutron stars at birth, slow (S) and fast (F) rotators. These 
two populations differ in their kinematical properties, and the F rotators are 
characterized by a mean velocity and by a scale height lower than the S ones. 
Basing on this picture, the time evolution of the distribution function in the 
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galactic potential has been studied by several authors (Blaes & Rajagopal 1991; 
BM; see also Paczynski 1990). In a recent paper, basing on a sample of 29 
young pulsars, Lyne & Lorimer (1994) suggested the possibility that neutron 
stars are born with typical velocities higher than both the F and S populations 
of Narayan & Ostriker. Their sample, however, is not complete and this result 
would prima facie imply that most pulsars evaporate from globular clusters and 
the Galactic plane. For this reason we retain the velocity distribution at birth 
proposed by Narayan & Ostriker. In any case, results presented by BM show 
that the contribution of high-velocity objects is really negligible as far as the 
observability of ONSs is concerned. 

Detailed tabulations of the distribution function are not available in the 
literature, so we have repeated the calculations presented by BM, evolving 
spatial and velocity distributions corresponding to the F population of model 
b of Narayan and Ostriker (1990). All NSs were assumed to be born at t = 0. 
The Galactic potential is sum of three contributions, disk, spheroid and halo, 
and is taken from Blaes & Rajagopal (1991). The fraction of F objects is 
55 % of the total number of neutron stars, and we use N tot = 10 9 . The 
distribution function was calculated integrating the orbits of 48000 stars, up 
to t = 10 10 ys. Initial conditions at t = were obtained assuming that the 
velocity distribution is gaussian in each component, together with the vertical 
distribution; the radial one is poissonian and the azimuthal one is uniform. 
Orbits were accepted only if the final parameters were in the ranges < R < 20 
kpc, \z\ < 2.5 kpc and \v\ < 200 km/s. We have found more convenient for later 
applications (see next sections) to store the distribution function in heliocentric 
coordinates, / = f(r,b,l,v). The intervals in heliocentric coordinates and in 
v were divided in 75 equally spaced bins. In order to compare the output 
of our simulation with the results of BM, we performed the evolution of the 
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distribution function for 20000 stars using galactocentric coordinates: in the 
local region 7.5 kpc < R < 9.5 kpc we find (z 2 ) 1 / 2 = 739 pc; (z) = 261 and 249 
pc in the northern and southern hemisphere, respectively. The mean velocity 
in this local region, averaged over \z\ < 200 pc, turns out to be 78 km/s, and 
(ij 2 ) 1 / 2 = 87 km/s. All these values are in close agreement with the results of 
BM. 

In the next sections, we will use /(r, 6, /, v) to estimate the observability of 
ONSs as X-ray sources. However, mainly in the case of spherical emission, a 
detectable X-ray flux at earth is produced only by the small number of stars 
that populate the low-velocity, small-r tail of the distribution. In this case a 
reliable statistical sample could be obtained only from the evolution of a large 
number of stars ( ~1% of the total population), requiring an unacceptably 
high computational time. In the alternative hypothesis in which the birth-rate 
of NSs is assumed to be constant, the statistics could be improved without 
increasing the number of computed orbits, as proposed by BM. In this case, 
in fact, it is convenient to store the distribution function every 10 6 ys: the 
final distribution can be obtained as the superposition of different outputs and 
will be representative of stars born at different times with the same initial 
conditions. In our case, however, we need to know the dependence of / on all 
spatial coordinates, and the storage and handling of such a high number, ~ 10 4 , 
of very large matrices, 75 x 75 x 75 x 75, become troublesome. 

To circumvent this problem, we have calculated the local velocity 
distribution g(v), averaged over the region r < 2 kpc, and a best fit Gfi t (v) of 
its cumulative function G(v) = J Q " g(v') dv'. We found that a function of the 
form 

GMV) = 1 + (v/v T (7) 
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with v = 69.0 km s -1 , n ~ m = 3.3, provides a good fit to G(v). Figure 2 
shows G(v) at t = 0, after t = 10 10 ys and Gf it . In this local region, the 
number density of stars turns out to be no = 3 x 10~ 4 (A^ O i/10 9 ) pc -3 , which is 
about 2.5 times lower than the midplane density found by BM. This apparent 
discrepancy is due to the fact that we are averaging over a sphere of radius 2 kpc, 
centered on the sun. By comparison the star density computed in BM at z = 1 
kpc is about 1 x 10- 4 (^ ot /10 9 ) pc" 3 . It can be seen that the relative number 
of low velocity stars is higher than the initial one, reflecting the fact that they 
do not move significantly away from the Sun during their evolution, as already 
noted by BM. The results of the numerical simulations are also confirmed by 
an alternative and faster semi-analytical approach based on an approximation 
of the 3th integral of the motion for low velocity stars (BM). 



4. THE CONTRIBUTION TO THE SOFT XRB 

The possibility that a new, yet undetected, population of X-ray sources 
contributes to the X-ray background, soothing the problem of the observed 
soft excess, has been recently suggested by Hasinger et al. (1993). By carrying 
out a detailed analysis of 27 fields at high galactic latitude, \b\ > 30°, these 
authors found that, at the faintest flux limit of ROSAT, about 60 % of the 
background is resolved into extragalactic discrete sources (see also Comastri et 
al. 1995). Their flux is 1.48 x 10" 8 erg cm" 2 s" 1 sr" 1 in the 0.5-2 keV band over 
a total flux Fxrb = 2.47 x 10 _8 erg cm~ 2 s _1 sr _1 in the same energy band. The 
projected number density of resolved sources turns out to be 413 deg -2 , which 
exceeds by about 60 % the density of QSOs predicted by standard evolutionary 
models. Even assuming that at the same flux limit the contribution of stars is 
~ 10 %, about 120 sources deg -2 remain to be explained, together with ~ 40% 
of the XRB (soft excess). As noted by Hasinger et al., these results can be 
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interpreted either by a more complicated model for the evolution of the QSOs 
X-ray luminosity function or invoking the presence of a new population of 
sources. An important point emerging from ROSAT observations that sets 
constraints on the nature of the new population is that the average spectrum 
of resolved sources becomes harder at lower flux limits. 

Recently MG discussed the possibility that the required new population is 
galactic in origin and derived its main properties, assuming that it is composed 
of standard candles of fixed luminosity L. In particular, they concluded that 
objects with a typical X-ray luminosity ~ 10 30 erg/s, and with a local density of 
~ 10~ 4 — 10~ 5 pc~ 3 could explain a substantial fraction, 20-40 %, of the XRB, 
giving, at the same time, the required source density of ~ 120 deg -2 . Assuming 
a blackbody spectrum at the star effective temperature, MG have ruled out 
the possibility that accreting ONSs could be the proposed population, just on 
the basis of their spectral properties. The emitted spectrum, in fact, is too 
soft, with a mean photon energy ~ 30 eV for L ~ 10 30 — 10 31 erg/s. The role 
of the magnetic field which drastically diminishes the emitting area, and the 
harder spectra obtained from our calculations, together with their hardening 
with decreasing luminosity, lead us to reconsider the possibility that ONSs can 
contribute to the XRB. 

Let us first consider the contribution of the entire ONSs population to the 
XRB. Here and in the following, we use the computed distribution function and 
the ISM density discussed in section 2 to evaluate the actual ONSs luminosity. 
The presence of the Local Bubble of radius ~ 100 pc surrounding the Sun 
has been neglected, since the fraction of stars that falls within this underdense 
region turns out to be negligible. The emitted flux per unit solid angle is: 
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where the integrals in r, O span the entire Galaxy and L x is the source 
luminosity in the 0.5-2 kev band, corrected for interstellar absorption. The 
upper limit in the velocity integral is set by the fact that no synthetic spectra 
are available for £ < 4 x 10~ 8 which corresponds to a minimum flux of 
3.8 x 10 17 erg cm _2 s _1 at the source. This is of little importance, since equation 
(5) shows that £ x becomes vanishingly small for l^2x 10~ 8 . 

Next we consider the number of resolved sources which is given by 



produces a count rate above the threshold of Hasinger et al. (estimated 
to be ~ 2 x 10~ 4 counts s -1 , corresponding approximately to a flux of 2 x 
10~ 15 ergcm _2 s _1 ). Both calculations have been repeated in the case of polar 
cap accretion with B = 10 9 G; the total column density corresponding to a given 
star position was used. The spatial distribution of ONSs has a characteristic 
scale height of ~ 260 pc, so both contributions, to the XRB and to the 
number count of sources, depend on direction. In order to quantify the degree 
of anisotropy, the two calculations have been repeated for objects with high 
(|6| > 30°) and low (|6| < 30°) galactic latitude. The first case allows a closer 
comparison with Hasinger et al. and MG. 

Our results are summarized in table 1. In the case of spherical emission 
the ONSs contribution to the XRB is really unimportant (< 1 %) similarly 
to the finding of MG & BM. Only a small fraction of stars, those in the low 
velocity tail of the distribution, have spectra hard enough to give a non negligible 
contribution to the background intensity. Our result is, anyway, ~ 40% greater 
than the estimate by BM. 

In contrast, when a non-zero magnetic field is taken into account, the ONSs 





where d. 



max 



is the maximum distance at which a star of luminosity i 
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contribution to the XRB, averaged over all latitudes, is found to be of order 
10%, and 5-6 % at the high latitudes considered by Hasinger et al. and MG. 
This corresponds to 12-25% of the observed soft excess. As can be seen from 
table 1, the ONSs contribution is higher by a factor 2-3 in the galactic plane, 
implying a certain degree of anisotropy. However, this is not in contrast with 
the observed near-isotropy of the XRB at energies in the 0.5-1 keV band (see 
e.g. Mc Cammon & Sanders 1990) since ONSs can contribute at most up to 
25% of the soft excess. At the sensitivity limit of the deep survey of Hasinger 
et a/., the number density of resolved sources, averaged over all latitudes, is 
~ 20 deg -2 , which amounts to about 5% of the total and to ~ 18 % of the 
non-QSOs, non-stellar component. 

From a direct analysis of the on line catalogue of the ROSAT PSPC 
pointings (ROSATSRC), we found that the mean number of sources detected 
with flux larger than 10~ 3 counts s -1 (~ 10~ 14 erg cm~ 2 s _1 ), is ~ 30 deg -2 
(b < 30°) and ~ 40 deg -2 (b > 30°). This result is formally in agreement with 
the mean number of unidentified sources predicted by our model at the same flux 
limit, ~ 10 deg -2 (see the next section), although it is apparent that the ratio of 
ONSs over the total number of sources is unacceptably high. The contribution 
of ONSs to the soft X-ray background would scale as the percentage of sources 
which will be identified as ONSs. 

The estimated ONSs contribution to the diffused emission is not probably 
strongly affected by the limits of the statistical analysis and the details of the 
ISM distribution, as it is shown below. Assuming a mean value n# = 0.5 cm -3 , 
B = 10 9 G and using equation (5), the X-ray luminosity is: 



In order to estimate the contribution provided by the most luminous stars, we 




(3.5 x 10" 4 - 5.5 x 10 _5 lni;) /v 3 , v>c 



s ? 
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can expand the approximate cumulative function, discussed in section 3, for 
(v/v ) n = x n C 1, to get (m ~ n) 

J v ' dv v 

The mean source luminosity is then 

P^max 

(4) ^ / £ x (v)f(v)dv, 
Jo 

where v max is the largest velocity at which our expansion can be reliably used. 
If we assume that the largest emission comes from stars within a volume V of 
radius 1 kpc with density n = 3 x 10~ 4 pc -3 and at a mean distance d ~ 500 
pc, the contribution to the XRB turns out to be: 

{(■x) n jr i no w ln -9 2 -1 „ -1 



-V = 1.23 x 10" y ergcm^ s" 1 sr" 1 ~ 0.05F X rb , 



And 2 4n 

that is of the same order of that obtained from the numerical simulation. 

5. OBSERVABILITY OF INDIVIDUAL SOURCES 

In this section we rediscuss the observability of accreting ONSs, following 
previous studies on this subject (TC; BM; CCT). In particular we refer to 
PSPC on board of ROSAT; the response curves, needed in equation (4) to 
obtain the count rate, are taken from ROSAT guide for observers. Due to the 
combined differential effects introduced by the response of the detector and the 
absorption of the ISM, this analysis crucially relies both on the spectral shape 
and on the mean photon energy. 

As can be seen from equation (2), typical values for n# are about 0.2 — 
0.5 cm -3 . Higher values seem to occur in the Sco-Cen and Per directions, 
where n# ~ 1 — 10 cm -3 . As far as the observability of individual sources is 
concerned, we assume a uniform medium and treat nu as a free parameter. This 
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will also allow a closer comparison of our results with those obtained in previous 
investigations (TC; Madau & Blaes 1994). We consider here three typical values 
for the density, tir = 0.2, 0.5, 1 cm -3 , as representative of different lines of sight. 

A star of luminosity L (i.e. moving at v = v) can be observed up to a 
maximum distance d max at which the count rate N (eq.[4]) becomes lower than 
the sensitivity limit of the detector (1.5 x 10~ 2 counts/s for the all sky survey, 
ASS, and 10~ 3 counts/s for Deep Exposure, DE). Table 2 shows d max for two 
typical luminosities, £ = 10~ 7 and £ = 4 x 10~ 8 : at distances d < d max all stars 
with L > L (i.e. v < v) can be observed. For comparison the corresponding 
quantities, as derived for a blackbody spectrum are listed in parenthesis. 

The actual distribution function f(r,v,l,b) of ONSs was derived in section 
3; since the scale height of the distribution, ~ 250 pc, is comparable to that 
of the gas, / can reasonably provide an estimate also of the actual distribution 
of accreting ONSs, so we could, in principle, calculate the expected number of 
observable objects within d max : 

N ons (<d max ,<v)= / / fr 2 drdQdv (10) 

Jn Jo Jo 

where r is the radial distance from the Sun. To avoid possible fluctuations due 
to the small number of slow stars in the local region present in our simulation, 
we replaced the previous formula with 

{^dmaxnoG fit (v) d max < zism, 

(11) 
max > zism, 

where n = 3 x 10~ 4 pc~ 3 . The typical scale height for the ISM is zjsm = 
300 pc; ONSs with z larger than zjsm are assumed not to accrete. The 
corresponding values are reported in table 2. Eq. (11) contains, nevertheless, 
two major simplifications: first, we assumed that locally the spatial and 
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velocity dependence of the distribution function can be factorized; second, 
we considered uniformly distributed ONSs. Since this could invalidate our 
statistical analysis, we stress that different columns in the table do not have 
the same reliability. A comparison with the corresponding values calculated 
for a blackbody shows that maximum distances computed for actual spectral 
distributions are systematically higher and in high-density regions can become 
greater by a factor 2-3. All calculations have been repeated for polar cap 
accretion with B = 10 9 G and for the same values of £ (see table 3). 

We note that stars with £ 10~ 8 are detectable up to typical distances 
~ 200 — 300 pc for ASS, even if they emit from the entire surface. A comparison 
with TC results shows that these values are of the order of their estimates for 
polar cap accretion. However, we stress that in our case the spectral hardening 
is only due to radiative processes; when a non zero magnetic field is taken into 
account, objects with the same luminosity become visible at larger distances. 
As it is apparent from table 3, if the magnetic field is 10 9 G about 10 sources 
deg -2 are expected to be above the DE ROSAT threshold; this figure is a factor 
~ 10 greater than the estimates of TC and CCT. 

It has been already stressed (BM; CCT) that giant molecular clouds in 
our galaxy provide the most favorable environment for the occurrence of high 
accretion rates. We reconsider this possibility, referring to the sample of 18 
nearby clouds studied by Dame et al. (1987). Assuming that the clouds 
are spherical and homogeneous, the values of their radius R c and density n c 
has been calculated by CCT. A luminosity £ = 10~ 7 , which correspond to 
typical velocities in the range 40 — 90 km s _1 has been assumed. Since velocities 
are largely supersonic, the use of equation (1) is fully justified. Being the 
local density much higher than the typical ISM density, we assume that all 
absorption is due to the cloud material. Table 4 shows the count rate (CR) for 

18 



the calculated spectrum and the expected number of observable ONSs: 

N ons (<v) = ^R 3 c n G(v); (12) 

the blackbody count rate is also reported for comparison. The count rate turns 
out to be above the threshold of Rosat (in DE) for more than half of the clouds 
in the spherical case, due to the hardening af the spectra with respect to a 
blackbody one. Clouds 15-16 and 18 appear also promising because of their 
high expected count rates, although the expected number of ONSs is small. 
We note, however, that, since the number of neutron stars in the clouds are 
small, large statistical fluctuations are possible. The corresponding values for 
polar cap accretion are summarized in table 5. Clouds 5-7, 9-13 represents the 
most favorable sites for observability and this agrees with the previous estimate 
presented by CCT. 

6. DISCUSSIONS AND CONCLUSIONS 

By computing the present distribution function of ONSs and using a map of 
neutral and molecular gas in the local region, we performed an analysis of the 
ONSs contribution to the diffuse emission in the soft X-ray energy band. A 
recent calculation of the spectrum emitted by accreting neutron stars (Zampieri 
et al. 1995) shows that, for luminosities ~ 10 30 — 10 31 ergs -1 , a relevant fraction 
(10%-20%) of the total flux is emitted in the 0.5-2 keV band even in the 
absence of magnetic field. Furthermore, the spatial density of ONSs in the solar 
neighbourhood computed using the evolved distribution function turns out to 
be no = 3 x 10~ 4 (iVt o t/10 9 ) pc~ 3 . Both these values are in agreement with the 
constraints set by MG on possible Galactic candidates which may contribute 
substantially to the XRB. We found that, in the case of polar cap accretion 
(B ~ 10 9 G), when a greater number of objects produce spectra hard enough 
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to be detectable, the ONSs contribution is « 10% of the total measured intensity 
in the 0.5-2 keV band, corresponding to 25-50% of the observed soft-excess. 
The fraction of resolved sources at the limit of the deepest ROSAT surveys is 
about 5 %. However, the strongest constraint follows from the expected number 
of sources above 10~ 3 counts s _1 and this will determine the importance of the 
ONSs contribution to the X-ray background. 

We note that these results are not strongly dependent on the overall shape 
of the velocity distribution of ONSs, but only on the number of objects in the low 
velocity tail, here calculated using the Narayan & Ostriker distribution function. 
For a high velocity distribution at birth, such as that recently presented by Lyne 
& Lorimer (1994), the extrapolation to low velocities seems to indicate that the 
number of stars with velocity less than ~ 30 km/s, turns out of the same order. 
However one should keep in mind that the distribution in the low velocity tail 
of ONSs remain very uncertain; for instance, it can be affected by dynamical 
heating, as suggested by Madau & Blaes (1994). This process, observed in the 
local disk star population, causes the velocity dispersion to increase with age as 
a consequence of scattering by molecular clouds and spiral arms (Wielen 1977). 
If ONSs participate the same process, dynamical heating over the lifetime of 
the Galaxy may scatter a fraction of low velocity stars to higher speeds. This 
could have a major effect on the source number counts (that can be decreased 
up to a factor 10) and it may reduce the contribution of luminous ONSs to the 
background. Other factors of indetermination, like the poor knowledge about 
the NSs birth rate, and hence, of their present total number, may affect our 
conclusions to the same extent. In addition, present estimates of the magnetic 
field of ONSs are subject by various uncertainties. 

We conclude that, even if neutron stars do not account completely for the 
characteristics of the galactic population proposed by Hasinger et al. (1993) 
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and MG, their contribution may be of importance. 
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FIGURE CAPTIONS 

Figure 1. Emergent spectra for L = 10~ 7 and 4 x 10~ 8 LEdd (full lines), 
together with the corresponding blackbody spectra at the neutron star 
effective temperature (dashed lines). 

Figure 2. G(v) at t = (dashed line) and after t = 10 10 ys (continuous line); 
Gfu (dash-dotted line) is also shown. 
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